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Abstract 

We investigate the effect of the deuteron breakup on the optical 
potential for the elastic scattering of 56 MeV deuterons from 51 V. The 
breakup probabilities calculated within the post-form distorted wave 
Born-approximation theory of breakup reactions are fitted to derive 
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the contribution to the optical potential from the breakup channels. 
The breakup potentials are found to have large real as well as imag- 
inary parts. Thus the dynamical polarization potential due to the 
breakup process is expected to modify strongly the real part of the 
optical potential calculated by the double folding model. 
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1 Introduction 



The folding model with realistic effective nucleon-nucleon interactions has 
given a good insight into the nucleus-nucleus interaction (see e.g. a recent 
review by Brandan and Satchler However, when a light ion interacts 
with other ions there is also a possibility of the breakup of the projectile into 
two or more fragments. If the breakup channel is strong, it will affect not 
only the imaginary potential but also the real one. This leads to a dynamical 
polarisation potential (DPP) which has to be added to the real potential 
obtained by the folding model. The DPP due to the breakup process can, 
for example, be estimated from the adiabatic model proposed by Johnson and 
Soper pi, or by more sophisticated continuum discretised coupled channels 
techniques M. The DPP is required to describe the elastic scattering of 
the weakly bound projectiles like 6,7 Li and 9 Be in the folding model with 
reasonable renormalisation factors 0]. 

Experimental studies have shown that the breakup probabilities increase 
drastically with energy even for tightly bound projectiles |4], |J. For example, 
the cross section for the breakup of the a particle increases by, at least, 
an order of magnitude as the beam energy is varied from 65 MeV to 140 
MeV H 0j. At higher beam energies (> 140 MeV) the breakup cross section 
can be as large as 25% of the total reaction cross section. Thus the DPP 
could be required also for tightly bound projectiles for beam energies above 
30 MeV/A. 
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In a recent report || one of us investigated the contributions to the a- 
particle optical potential from the (a, 3 He) breakup reaction on 62 Ni target 
at the incident energy of 172.5 MeV. The calculations proceed in two steps. 
First the breakup probabilities are calculated within a theory which is formu- 
lated in the frame-work of the post form distorted-wave Born-approximation 
(PFDWBA). This theory has been found to reproduce the experimental data 
on the breakup of light projectiles extremely well || . In the second step, these 
probabilities are fitted to generate the breakup part of the optical potential 
which is assumed to consist of two parts; one which is entirely due to the 
breakup of the projectile and the another independent of it. It was found in 
[§] that breakup contributes substantially to the optical potential. However, 
in ||, the DPP for the alpha elastic scattering was obtained by fitting only 
to the (a, 3 He) breakup channel. For a complete determination of the DPP 
due to the breakup process the breakup probabilities for all the a-breakup 
channels, namely, (a,p), (ot,n), (a,d), (a, 3 He) and (a,t) should be calculated 
and fitted. 

Before taking up this rather ambitious task, which nobody has inves- 
tigated so far, we considered it worthwhile to perform calculations for the 
simple system of the deuteron to test the method in detail. In this case there 
are only two breakup channels, (d, p) and (d, n) and the experimental data on 
the deuteron breakup (see e.g. UlQ| ) is quite comprehensive. In this paper, we 
investigate the scattering of 56 MeV deuterons from 51 V. The breakup prob- 
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abilities for the (d, p) and (d, n) channels are calculated for this particular 
case. These are fitted to determine the DPP due to the breakup process. 

In sections 2 we describe the calculation of the breakup probabilities in 
the optical model. The calculation of the same within the PFDWBA theory 
of breakup is presented in section 3. The results for the DPP obtained by 
fitting the total and elastic breakup probabilities are discussed in sections 4. 
Our conclusions are presented in section 5. 

2 Calculations of breakup probabilities from 
the optical model. 

The calculations are based on the same principles as in Ref. ||. We assume 
that the optical potential V(r) is known from the phenomenological analyses 
of elastic scattering data and solve the radial Schrodinger equation 

_J^L± + [ k 2 _ [/( r ) _ w ( r ) = o, (l) 

where k is the wavenumber and U (r) is related to the optical potential V(r) 
through U(r) = (2m/h 2 )V(r). The solutions for the radial wavefunctions, 
yt{r), are normalised according to 

yt(r) = e i5e [cos 5 e F e (kr) + sm6 e G e (kr)}, (2) 

where Fi and Gi are the regular and irregular Coulomb functions. With this 
normalisation the partial wave amplitudes, T#, are written in terms of the 
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phase-shifts 5e as 

T e = e iSe sm5 e (3) 
The partial wave amplitudes may also be calculated from the relation 

T t = -k~ l / F e (kr)U(r)y e (r)dr (4) 
J o 

Now we define the optical potential due to breakup as Ub u (r) and write 
the full potential U(r) as (U(r) — U^r)) + Uf, u (r). With this definition, Eq. 
1 can be recast as 

+ [k 2 - (U(r) - UUr)) ~ U^r) - ^J^] y e (r) = (5) 
A discussion of this problem may be found in Ref. [11], [12] . 



By solving the Schrodinger equation for the "bare" potential U(r) — 
U bu (r), 

f^l + r . mr) .U^.^ ]Mr ^ , (6 ) 

we obtain the radial wavefunctions vi(r), the phase shifts S' e , and the partial 
wave amplitudes, T[. 

The relation between the partial wave amplitudes are 

T t = T[ + T\ u (7) 

The partial wave amplitudes for breakup are thus given by the difference 
between the partial wave amplitudes obtained for the potentials U(r) and 
U(r) — Ubu(r), respectively. These can also be calculated from the formula 
(see e.g. 0, |TJ) 

6 



roo 

T b t u = -AT 1 / v e (kr)U bu (r)ye(r)dr (8) 
Jo 

It should be noted that the integral in Eq. 8 contains the radial wave- 
functions ve, obtained with the bare potential U(r) — Ub u (r), as well as the 
radial wavefunctions ye, obtained with the full potential U(r). The assump- 
tion implicite therein is that the bare potential cannot lead to the breakup 
reaction. 

The breakup probability for a certain £-value is determined from the 
partial wave amplitudes Tf 1 according to 

Breakup probability = |T^ U | 2 (9) 

These breakup probabilities are compared with those calculated within 
the PFDWBA theory of breakup reactions (discussed in the next section) to 
determine the breakup potential, Ub u . In our procedure, this potential (with 
certain apriori assumed form) is varied to reproduce the breakup probabil- 
ities calculated within PFDWBA. It may be remarked here that when the 
potential Uf, u (r) is varied, new values of T[ have to be calculated from the 
solution of eq. 6, whereas the values of remain unchanged. 

It is interesting to note that the reaction cross section in the elastic scat- 
tering is known from the solutions of eq. (1). Furthermore the total breakup 
cross section can be calculated from the breakup probabilities. Even then, 
the reaction cross section for the bare potential U(r) — £/ 6u (r), will depend 
on the shape of the potential, £4 u (r). 
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In Ref. ||, the peripheral dominance of the breakup process made us 
parametrize the breakup potential as a Gaussian with the addition of a 
Woods-Saxon form factor. In this work, however, this preconceived view 
is avoided and a Fourier-Bessel (FB) expansion 

U bu (r) = E a n W^, (10) 

n=l H 

is used for both the real as well as imaginary parts of the breakup potential. 
The radius R was kept fixed at 18 fm in all calculations. 

3 Microscopic calculations of breakup prob- 
abilities. 

In the theory of the breakup reaction (d — > p + n ) formulated within the 
framework of the post form distorted-wave Born-approximation, the proba- 
bility of breakup P\ up<yd,pS} is defined by IT 
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vtotal( d >P) = J dflpdEp 



dVlpdEp 



where ^f^g^ is the double differential cross section for the inclusive breakup 
reaction (d,p), which is the sum of the elastic and inelastic breakup modes. 
These are given by 

d 2 o (elastic) HdHpUn kpkn . rp 12 
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where the T-matrix Te 



Ti nmn = ^o/^ X (-)*(k p ,^r)^^ m „(f) X i +) (k d ,r)A(r)P(rj(13) 

-Do is the well known zero range constant for the d — > p + n vertex. A(r) 
and P(r) are the finite range and nonlocality correction factors respectively. 

are the optical model wave functions in the respective channels with k's 
being the corresponding wave vectors. 

The inelastic breakup cross section is given by 

d a [inelastic) _ ^-d^p^n kp^n sr^ (reaction /-elastic} I rp _ T° (&A\ 

dn p dE p - (2tt)^ 6 k d " 1 ' 1 nmn tnmJF ' ' 

In this equation (j^ actlon an d af^ stlc are the reaction and elastic scattering 
cross sections for the neutron - target system corresponding to the partial 
wave i n respectively. The T matrix T® nTnn is defined in the same way as 
Eq. (13) with the elastics scattering wave function xe n being replaced by the 
spherical Bessel function. It may be noted that Eq. (14) includes contribu- 
tions from all the inelastic channels of the neutron + target system ( see e.g 
H for complete detail ). 

The angular integration in Eq. (13) is performed by introducing the 
partial wave expansion of the distorted waves and using the orthogonality of 
the spherical harmonics. The resulting slowly converging radial integrals are 
evaluated very effectively by following the contour integration technique of 
Vincent and Fortune JTj| . 

We require the optical potentials in the incident and outgoing channels as 
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input in our calculations. In the results presented in this paper, the deuteron 
optical potentials were taken from the global sets given by Daehnick, Childs 
and Vrcelj fl4| whereas the potentials of Becchetti and Greenlees ]15j were 
used in the neutron and proton channels. 

The total breakup probability Pe d is defined as following: 

P b t; upid ' n) (inelastic) (15) 

In Fig. |l] we show the results for the breakup probability for the deuteron 
incident on a 51 V target at the beam energy of 56 MeV. We can see that 
the (d, p) and (d, n) breakup probabilities are similar in shape and absolute 
magnitude. The elastic breakup probability is much smaller and shows a 
different behaviour function of i*. 

As discussed in Ref. [[5| the total cross section can also be expressed as 

a breakup = 27T J db b P^ 4 (16) 

where the impact parameter, b, is related to the angular momentum, id, and 
the wavenumber kd, through b = (id + 1/2) 

The cross sections for the inelastic (d,p), inelastic (d,n) and elastic (d,pn) 
breakup processes were found to be 290 mb, 214 mb, and 122 mb respectively. 
This leads to a total breakup cross section of 627 mb. It may be noted that 
our total (d,p) breakup cross section (which is the sum of inelastic (d,p) and 
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elastic (d,pn) cross sections) is 412 mb which is in in reasonable agreement 
with the measured value of 481 mb reported by Matsuoka et al. [ID]. 



4 Results and discussions 

The fitting procedure was the same way as that described in Ref . |§ . How- 
ever, the imaginary part of the breakup potential, Ub u , was assumed to be less 
than that of the full potential U (r), so that (U (r) — f/ bM (r)) was absorptive for 
all the radii. The errors associated with the calculated breakup probabilities 
in the optical model (Eq. (9) ) were about 10%. In the fitting procedure, 
14 coefficients ( a n as defined in Eq. 10) were varied simultaneously for each 
potential to get a minimum in the \ 2 . 

In the calculations for (a— 3 He) breakup || it was found that the gross 
properties of the breakup probabilities could be reproduced with a purely 
real as well as a purely imaginary breakup potential. An imaginary part of 
the potential was needed for fitting them at only small f-values. However, 
in case of the deuteron we found it impossible to reproduce the breakup 
probabilities (Pe d ) without a complex breakup potential. One of the reasons 
for this could be the fact that the breakup probabilities in the this case are 
considerably larger than those for the (a— 3 He) reaction. Fig. [| shows the 
best fit obtained with a purely real breakup potential whereas Fig. ^ with 
a purely imaginary one. As can be seen from these figures, the imaginary 
potential is very important for reproducing the breakup probabilities at small 
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^-values, while the real part is required to reproduce these for the grazing 
partial waves and beyond. This suppports the observation made in [I7| that 
the real potential, due to non eikonal effects, increases the absorption at large 
radii and enhances strongly the peripheral collisions. 

The best fit result obtained with a complex potential is shown in Fig. f|. 
Both potentials include 14 terms in the FB expansion. It can be seen that the 
breakup probabilities calculated within PFDWBA (shown in Fig. |l|), are re- 
produced very well. The imaginary part of the breakup potential is found to 
be very strong. This is not surprising since the inelastic interactions (which 
could cause the breakup of the projectile) are taken into account by the imag- 
inary part of U(r). These will automatically be included in the imaginary 
part of Uf, u . In the folding model calculations of the elastic scattering where 
the imaginary part of the potential is phenomenologically determined, it is 
mainly the real potential which is of interest. Therefore, our result suggest 
that the inclusion of the potential, Uf, u , in these calculations will lead to a 
strong modification to the real part of the folding model potentials. 



Yabana et al |L8| have investigated the effects of breakup for the elastic 
scattering of 80 MeV deuterons from 58 Ni. Both real and imaginary parts 
of the DPP obtained by these authors are weaker than those obtained in 
our work. It may be remarked here that, these authors have neglected the 
Coulomb breakup process, which can be quite large for the deuteron target 
interaction [F20|, 01. In our work, on the other hand, this is included on the 
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same footing as the nuclear breakup. Moreover, our breakup probabilities 
lead to the cross sections which agree well with the experimental data on the 
deuteron breakup. 

We should, however, stress that, there are (as in optical model calcula- 
tions) ambiguities in the breakup potentials obtained by us. For instance, 
in Fig. H we show the results obtained when the number of terms in the 
real potential is decreased from 14 to 4. The fit to the PFDWBA breakup 
probabilities is still reasonable. This therefore, makes it difficult to arrive 
at a definite conclusion about the shapes of the DPP. Nevertheless, it is in- 
teresting to note that the real part of the DPP obtained in this way is still 
larger than those of Yabana et al. |18| 

Next we discuss the breakup probability for the elastic breakup ( a process 
in which the target nucleus remains in the ground state). In theories like the 
one suggested by Bertsch, Brown and Sagawa |jn| the reaction cross section 
in nucleus-nucleus collisions are calculated from collisions between nucleon- 
nucleon pairs in the projectile and target nuclei. In such approaches it is 
necessary to include a correction due to the elastic breakup of the projectile. 
We term the corresponding potential as the "dissociation" potential {Udis) to 
separate it from the potential Ub u defined earlier. It would also be worthwhile 
to study how the dissociation affects the optical potential and the reaction 
cross section. Udis is obtained by fitting (in the same way as described above) 
the elastic breakup probabilites as calculated by PFDWBA theory (see Fig. 
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|1|). We stress that these include both Coulomb and nuclear breakup as well 
as their interference terms. 

The results obtained with a purely real U^is are shown in Fig. || The 
fits to the elastic breakup probabilities are satifactory for the grazing partial 
waves. However, those for the lower partial waves are poorer. On the other 
hand attempts to fit them with a purely imaginary Udis resulted in a very 
bad agreement. Therefore we reduced the real potential shown in Fig. |6| 
and tried to reproduce the data by a variation of the imaginary potential. 
The results from this search is shown in Fig. [7| and it is evident that 
the maximum is very badly reproduced. However, when both the real and 
imaginary potentials are varied a good fit to the elastic breakup probabilities 
are obtained, even if the shape of the imaginary part of Ub r so obtained looks 
somewhat unusual. Figure |8| shows one of the best fit potentials obtained in 
this way. 

As discussed in section 2, the partial wave amplitude for the potential 
U(r) is given by the sum of the partial wave amplitudes the potentials, U(r) — 
Ub u (r) and Ub u (r). This is not true for the reaction cross sections calculated 
from each potential separately, since these are sensitive to interference effects. 
The reaction cross section for the full potential is 1571 mb. With the real 
potential shown in Fig. ^ the reaction cross section for the potential, U (r) — 
Ubu{r), is found to be 1405 mb. The difference (166 mb) is somewhat larger 
than the value (122 mb) obtained for the total elastic breakup cross section 
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using Eq. (11). With the complex potential shown in Fig. [8], the reaction 
cross section without breakup potential does not decrease, instead it goes to 
a value of 2173 mb. These examples show the importance of including the 
effects of the dissassociation in the optical model. The inteference effects are 
important and have to be treated correctly. 

Our calculations, therefore, indicate that the dissassociation of the deuteron 
should give a substantial contribution to the dynamical polarisation poten- 
tial. The assumption that the optical potential, obtained from a fit to elastic 
scattering data, include effects of dissassociation and inelastic breakup in 
a correct way must be questioned. There are examples, as in Refs. [|l"3|1 , 
where the absorption due to breakup is even larger than the total absorption 
predicted by the optical model for large ^-values. Since the elastic breakup 
is considerably larger than the inelastic breakup for large ^-values, this in- 
dicates that the optical models should be modified with a contribution of 
the dissassociation process, which may not necessarily have an effect on the 
angular distributions for elastic scattering. 

5 Conclusion 

In conclusion, our study of the effects of the breakup on the optical potential 
for the elastic scattering of 56 MeV deuterons from 51 V shows that breakup 
gives a substantial contribution to this potential. We found that the strong 
enhancement of the peripheral collisions requires that this contribution has 
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a real part. Thus the folding model calculations should include a dynamical 
polarisation potential. 

The investigations of the effects of the dissassociation also indicate a con- 
tribution to the dynamical polarisation potential. However, the conventional 
shapes used for the optical potential do not account for the dissassociation 
process in a correct way. Therefore the investigation should be repeated with 
more sophisticated optical model potentials with different shapes. 

We plan to calculate breakup probabilities for all breakup channels for 
a-particles. We believe that these together will give total breakup probalities 
of the same order as those of the deuterons studied here, at least for energies 
above 40 MeV/A. We see no reason why the breakup potential for a-particles 
should be purely imaginary when it is complex for deuterons. Therefore we 
believe that all light particles require dynamical polarisation potentials in 
folding model calculations at higher energies and that the energy depen- 
dence of effective interactions, which presently reproduce a-particle scatter- 
ing, should to be modified. 
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Figure 1: Calculated breakup probabilities in the scattering of 56 MeV 
deuterons from 51 V. 



Figure 2: Results obtained with a real breakup potential. In the upper 
part, (a), the dashed curve shows the nominal potential, the solid curve the 
breakup potential and the dotted curve the difference. In the lower part, (b), 
the open circles show the total breakup probabilities from Fig.l and the solid 
circles the fitted values 
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Figure 3: Results obtained with an imaginary breakup potential. In the 
upper part, (a), the dashed curve shows the nominal potential, the solid curve 
the breakup potential and the dotted curve the difference. In the lower part, 
(b), the open circles show the total breakup probabilities from Fig.l and the 
solid circles the fitted values 



Figure 4: Results obtained with a complex breakup potential. In the upper 
parts ,(a) and (b), the dashed curve shows the nominal potential, the solid 
curve the breakup potential and the dotted curve the difference. In the lower 
part, (c), the open circles show the total breakup probabilities from Fig.l 
and the solid circles the fitted values 



Figure 5: Results obtained when the number of terms in the FB-expansion 
of the real potential has been decreased from 14 to 4. The notations are the 
same as in Figure 4. 



Figure 6: Results obtained with a real dissociation potential. 



Figure 7: Results obtained with an imaginary dissociation potential. The 
real potential is in this case 50 % of the real potentiatial shown in Figure 6 



Figure 8: Results obtained with a complex dissociation potential. 
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